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Abstract 

Using a long wavelength iteration scheme to solve Einstein's equations near the Big- 
Bang singularity of a universe driven by a massive scalar field, we find how big initial 
quasi-isotropic inhomogeneities can be before they can prevent inflation to set in. 



PACS numbers: 98.80.Cq, 04.50.+h 



1 



1. Introduction 

The scale factor a(t) of a spatially flat isotropic and homogeneous Robertson- Walker 
(fRW) universe driven by a massive scalar field ip evolves as follows: near the Big-Bang 
(chosen to be t = 0) the scalar field goes to oo as — In t and behaves like a perfect stiff 
fluid (whose pressure equals the energy density, that is whose adiabatic index T is 2), so 
that a grows as t 1 / 3 [1]. Then a grows quasi exponentially during the inflationary regime 
during which ip slowly rolls down its potential well [2] . Finally at the end of inflation when 
ip oscillates in the bottom of the well a behaves as if the universe was driven by dust and 
grows on average as t 2 / 3 [3]. 

A question is : how stable is this evolution against departures from spatial flatness, 
isotropy, and homogeneity. The effect of curvature was studied in [4]; since its role af- 
ter the Big-Bang can become predominant, it can prevent inflation from ever starting if 
strong enough. Inflation in homogeneous albeit anisotropic Bianchi models was thoroughly 
analyzed, with similar conclusions : see e.g. [5] for a review. Finally, the role of inhomo- 
geneities was studied analytically under simplifying assumptions in [6] . They were studied 
numerically, in the case of planar symmetry in [7] and in the case of spherical symmetry 
in [8]. In [8], setting the initial conditions in such a way that inflation would occur only in 
the central region if the universe evolved like a fRW space time, the conclusion was that 
inflation would indeed occur only if the central region was larger than a few times the local 
Hubble radius. 

Here we shall tackle the problem semi-analytically in the long wavelength approxima- 
tion. 

The long wavelength iteration scheme, the history of which goes back to [9], is a way 
to build, out of "seed" spatial metrics, approximate solutions of Einstein's equations which 
describe inhomogeneous but quasi-isotropic universes on scales larger than the local Hubble 
radius (see [10] and references therein for a detailed description of the scheme). When 
matter satisfies the strong energy condition (that is does not inflate) this approximation is 
well suited to describe the early universe since most scales are then larger than the Hubble 
radius on average (this is the well-known horizon problem). These approximate solutions 
however are not generic in the sense that they are built out of a seed spatial metric, that is 
on three physically distinct arbitrary functions, instead of 6 (4 for the gravitational field, 
plus 2 for the scalar field) . The three missing functions can be seen as describing departures 
from isotropy. These anisotropies cannot be neglected near the Big-Bang (see e.g. [9-10]) 
but they decay with time much faster than all other contributions to inhomogeneity and 
hence will be ignored (see [11] for an analysis of the generic solution near the Big-Bang). 

Within that scheme the growth or decay of inhomogeneities according to their equation 
of state can be easily inferred, at least qualitatively (see [10, 12]): they decay when matter 
violates the strong energy condition, that is inflates, and grow otherwise, that is when the 
effective adiabatic index T e ff of matter is > 2/3. Therefore when matter is a scalar field 
the inhomogeneities first grow (T e ff = 2), then decay during the inflationary period, to 
grow again at the end of inflation (T e ff = 1), at the condition that they do not grow so 
large during the first phase as to prevent inflation to start. 

The purpose of this paper is to give quantitative estimates on when inflation may 
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be halted by the presence of initial quasi-isotropic inhomogeneities. This will be done by 
integrating numerically the ordinary second order linear differential equations that govern 
their evolution in the long wavelength approximation scheme. 



2. The equations 

The long wavelength approximation scheme [10] consists in looking for solutions of 
Einstein's equations for gravity coupled to a scalar field whose 3-metric (in the t = Const 
slices of a synchronous reference frame) can be expanded as a sum of spatial tensors of 
increasing order in the gradients of a "seed" metric with time dependent coefficients. The 
line element is thus of the form: 

ds 2 = ~dt 2 + 7y(t, X k )dx i dx j , =W 7y +^ 3) 7ij +^ 5) 7 y + ... 

^ ij = a 2 (t)h ij (x k ) , Wi ij = a 2 [a 2 (t)Kh ij + b 2 (t)R ij ] (1) 

(5) 7y =« 2 [a^t)R 2 hij + b A {t)RR l3 + c^t)R lm R lm h l3 + d A {t)R im Rf] 
+ a 2 [e±{t)V m V m RK 3 + f^V^.R + <7 4 (*)V m V m i2y] , 
and the scalar field of the form: 

<p = <p(t) + Mt)R + M*)R 2 + Mt)RimR lm + &(*)v m v m fl + ... (2) 

This solution is built on a "seed" hij(x h ) of 6 functions of space which can be reduced to 3 
by a suitable choice of spatial coordinates. As shown in [10], it is not a generic solution (it 
should then depend on 6 physically distinct functions) but it is an attractor of a class of 
generic solutions. At lowest order the metric reduces to 7^ = a 2 (t)hij(x k ) and cp = <f>(t), 
hence the name "quasi-isotropic" given to the solution. Vj is the covariant derivative with 
respect to hij, Rij its Ricci tensor, R = R^hij the curvature scalar. The coefficients 
a(t),a,2(t) etc are functions of time which are determined by Einstein's equations. We 
shall denote by L the characteristic comoving length on which the spatial metric varies: 
dajk ^ L~ 1 ^ jk . 

At zeroth order Einstein's equations reduce to the Friedmann equations for a spa- 
tially flat Robertson- Walker (fRW) universe and determine a(t) and <j>(t). Introducing the 
dimensionless variables: 

T = mt , F = 2v / 3^> , H = £ , S = am 1/3 (3) 

m a dt 

where m is the mass of the scalar field, they read: 

dS_ HS_ ^i = _ (H 2_ F 2 ) ^L = _ V / H 2_ F 2 (4) 

dT 3 ' dT K ' ' dT V • K } 

The solutions of eqs (4) depend on three integration "constants": the time T{x k ) of the 
Big-Bang that we shall restrict to be T = (see [11] for an analysis of delayed Big Bang 
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solutions); the size S(x k ) of the scale factor at some given initial time: a "constant" which 
can be absorbed without loss of generality in a redefinition of the seed metric h^; and 
finally the initial value F for the scalar field. A priori F depends on space but the first 
order Einstein equations (see [10] for details): diH = —4(dF / dT)diF imposes that it does 
not. Its value, the same for all x k , tells us which curve of the (F, dF/dT) phase diagramme 
the solution follows, and hence determines the total amount of inflation. We shall therefore 
integrate these equations up to the beginning of inflation with, as initial conditions at time 
T = e: 

e 1/s , H in = - , F in = F, (5) 



e being chosen such that (lne) 2 is numerically negligible compared to e -2 . If matter were 
a perfect fluid with adiabatic index T the scale factor S would grow as T 2 / 3r . In the 
initial regime (5) the scalar field hence behaves like a stiff fluid (T = 2) [1]. When F 
is large enough the solution then enters an inflationary phase characterized by a slow 
linear decrease of F and a quasi exponential growth of S. The comoving Hubble radius 
Lh = 3/ SH hence first increases as T 2 / 3 during the stiff fluid regime and then decreases 
exponentially. The moment it reaches its maximum value can be taken as the beginning 
of inflation. We shall denote L 1 ^ and F in f the values of the Hubble radius and the scalar 
field at that moment. 

At third order Einstein equations determine the time dependent coefficients a 2 (t), b 2 (t) 
and 02 (t) in (1-2). In terms of the variables defined in (3) and introducing: 



B 2 = b 2 m 4 / 3 , A 2 = a 2 m^ 3 , F 2 = A\j ^-m 4 / 3 <f> 2 (6) 



they can be written as (see [10]): 



du _ dB 2 _ u 

dr~ 2b ' It's* {7) 



dv 
dT 



u d 
4SdT 

Fo = 



S 



dF 
dT 



dA 2 
dT 



dF_ 

dT 



1 (ldB 2 dA 



dF/dT V 



4 dT 



+ 



dT 



(8) 



(9) 



It is easy to see that when the inflationary regime has set in, B 2 tends to a constant 
(that we shall call I 2 ) as well as F 2 . As for A 2 it decreases linearly in time. Now, as 
shown in [10], a gauge transformation modifies the coefficient of the scalar curvature R 
in ( 3 "*7ij but leaves untouched the Ricci term. During inflation it amounts to adding an 
arbitrary linear function of time to A 2 . Therefore the linear decrease of A 2 can be gauged 
away by a suitable choice of initial conditions for A 2 , v 2 , leaving B 2 as the only relevant 
quantity to be studied. Physically the geometry evolves from one configuration, hij, to 
another, + (I 2 /m 4 ^ 3 )Rij, where I 2 is the "imprint" on the geometry left by the initial 



4 



inhomogeneity. We shall therefore integrate (7) with the initial conditions (5) at T = e 
together with dF/dT\ in = — 1/e and: 

u in = -\e^ , B 2 \ in = - 9 -e*/\ (10) 



3. The solutions 



Equations (4) with initial conditions (5) and equations (7) with initial conditions (10) 
are integrated numerically. The results can be encapsulated in a plot of r 2 = a/I-T^/L^ 
as a function of Fi n f. See the curve k = of Fig 1. We see that as i^ n / increases, r 2 
tends to a constant close to 1 (when a brutal extrapolation of the analytical behaviours 
(5) and (10) would have given l/2y/2). Going back to the expressions (1-2) for the metric 
and the scalar field, a sufficient condition for inflation then appears clearly : inflation will 
set in if the corrective terms ^ '^iij remain small compared to ' '7^ that is if I 2 /L 2 < 1 
which is equivalent, because r 2 tends to 1, to L > L 1 ^ . This condition says that the size 
L of the inhomogeneity must be larger than the (comoving) Hubble radius at the onset of 
inflation, when the brutal analytic extrapolation gives that it can be smaller by a factor 
2\f2 = 2.8. Moreover, when the zeroth order barely inflates, that is when F in f is small, 
T2 = a > 1, so that a sufficient condition for inflation is L > aL 1 ^ . We thus recover in 
this semi-analytical approach the result of ref [8] . 



4. Convergence of the iteration scheme 



A rigourous mathematical analysis of the convergence of the series (1-2) is certainly 
beyond the scope of this paper. Using the results of the appendix of ref [10] we can however 
compute the next order, that is the coefficients 04,64 etc and see if including them spoils 
the conclusion of the preceeding paragraph, drawn from the first iteration. 

The relevant coefficients, which are not affected by a gauge transformation, are C4, 
c?4 and (74. We shall concentrate on g& which satisfies the following linear second order 
differential equations: 

dw dG 4 w 

df- B2S ' lT-'S^ 

where we have introduced G4 = m 8 ^ 3 g4 with the initial conditions: 

^ = -f/^ > G 4 \ in = - 8 ^. (12) 

(Knowing G4 and B 2: the coefficient D4 = rrft^d^ is readily obtained : D 4 = —\G 4 + \B\\ 
see [10].) 

The results of the integration of (11-12) are summarized in the diagramme r 4 = 
1 h as a function of F in f (see the curve k — of Fig. 2). We see that T4 follows 
the same pattern as r 2 . As F in f increases it tends a constant value of the order of 0.63 
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when an extrapolation of the analytical behaviour (16) would have given 8 3 / 4 = 0.21. As 
for the ratio | j D 4 | 1/4 /^h / it tends to 0.60. 

The inclusion of these fourth order terms therefore does not spoil the conclusion of the 
previous paragraph and, since T4 < T2, is an indication (although certainly not a proof!) 
that the exact solution does not differ drastically from the second order solution. 



5. Curvature versus gradient effects and an improved scheme in the case of 
spherical symmetry 



In the case when the inhomegeneity is imposed to be spherically symmetric, the seed 
line element da 2 = hijdx l dx^ can be written in the form: 



da 2 = 



dr z 



1 — kr 2 



+ r 2 dQ 



(13) 



where dfi is line element on the sphere and where k(r) is an arbitrary function of r. 

Expressing the Ricci tensor and its derivatives in function of k(r) we can write the 
long wavelength metric (1) and scalar field (2) under the form: 



ds 2 = -dt 2 + R 2 {C 2 dr 2 + r 2 dfl) 



with 



D _ r>app 
n — ^cRW 



I + j a 2 + ) k'r + 



Krw = a [1 + (3a 2 + b 2 ) k + 



C = C C RW 



1 



! h ( a 2 + -62 I k'r + ... 



, C c rw = 



and 



VT — kr 2 

P® + 2^ 2 + ... , ^ = + 6fcfc + ... 



(14) 

(15) 
(16) 
(17) 



where a prime denotes differentiation with respect to r and where the time dependent 
coefficients (a, 0), (a 2 , 62,^2) etc are the same as before and satisfy eqns (3-5), (6-10), 
(H-12). 

The rewriting of eq (1-2) under the form (14-17) shows clearly the two ways in which a 
spherically symmetric inhomogeneity makes the solution depart from the flat fRW solution. 
The first (trivial) effect is that of curvature: if k 7^ but all its derivative are taken to 
be zero, the long wavelength solution (14-17) can be checked to be nothing but the Taylor 
expansion in t of the exact curved Robertson- Walker (cRW) solution. In this case then 
the exact metric and scalar field are given by (14) with C = (1 — /cr 2 ) -1 / 2 , and R — a and 
4> satisfying the Friedmann equation: 



dS HS 



dH 



6K 



dF 



dT 



3 • rfT=- (if2 - F2 »-^ ' W —y/H*-*° + »K/S>. (18) 



where T, F, H, S are defined by (3) and where K = km 2 / 3 . 
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Integrating (18) numerically with the same initial conditions as before (eq(5)), we 
recover the results of ref [4], that is that if the curvature term K/S 2 is too large then 
inflation is halted. This can be seen in Fig. 3 where Lh is plotted as a function of T for 
different values of K and a given initial value F for the scalar field. We see that a negative 
K "favours" inflation whereas a positive value delays its setting in and, if large enough, 
can even prevent it (dotted line). 

The second effect is that of the gradients, that is the derivatives of k(r), which reflects 
the point to point correlation due to the variation of the curvature. It can be enhanced 
by choosing a seed k(r) such that it is small everywhere but has a steep gradient k'r 
around, say, r = R. From (16) and the result of section 3 that r 2 ~ 1 we know that these 
gradient effects will not prevent inflation to set in if the function k(r) is everywhere such 
that k'r < 2m 4 / 3 /(L~ / ) 2 . 

One can also improve the long wavelength scheme by replacing in (14-17) the ap- 
proximate cRW values by their exact values as given by (18) and taking for A2, -82,-^2 
the solutions of eq (7-9) where S,H and F are taken to satisfy (18) instead of (4). The 
results are summarized in Fig. 1 and 2 and show that the improved and the standard 
schemes coalesce for large Fi n f that is for strongly inflating solutions. This confirms what 
the previous section already indicated, that is that the iteration scheme seems to converge 
nicely. 

6. Conclusions 

An important question in inflationary cosmology is: how generic is it? As already 
shown by Goldwirth and Piran and confirmed here, inflation by itself requires a certain 
level of homogeneity: it can start if initial inhomogeneities are larger than the local Hubble 
radius. While the numerical calculations of [8] can explore in a detailed way a specific 
(spherically symmetric) case of strong initial inhomogeneity, the semi-analytical approach 
presented here is limited to rather small perturbations but it gives a better global picture 
on the factors that control the behaviour of the system. It gives only sufficient conditions 
for the onset of inflation, not as strong as the necessary conditions obtained in [8] but they 
are general and do not assume any spatial symmetry. We plan to extend the comparison 
between these analytical and numerical approaches by solving an identical initial value 
problem. This requires further investigation of the various coordinates systems used. 
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Figure Captions 

Figure 1 Plot of r 2 as a function of F in f. The full line corresponds to k = 0, the 
dashed line above to negative k, the dashed dotted lines below to positive k. 

Figure 2 Plot of T4 as a function of Fi n f. The conventions are the same as in Fig. 1. 

Figure 3 Plot of the comoving Hubble radius as a function of time for various values of 
the curvature term k and a given initial value for the scalar field. The full line corresponds 
to k = 0, the dashed lines below to negative k, the dashed dotted ones to positive k. 
Finally the dotted line is an example of a curvature strong enough to prevent inflation 
from setting in. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/9409056vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/9409056vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/9409056vl 
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